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We explore the different possibilities for branons as dark matter candidates. We consider a 
general brane-world model, parametrized by the number of extra dimensions TV, the fundamental 
scale of gravity Md, the brane tension scale / and the branon mass M. We analyse the parameter 
region in which branons behave as collisionless thermal relics (WIMPs), either cold or hot (warm) 
together with less standard scenarios in which they are strongly self interacting or produced 
non-thermally. 
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PACS numbers: 95.35.+d, ll.25.-w, ll.10.Kk 
INTRODUCTION 

Cold dark matter (CDM) is a fundamental ingredient 
of the current cosmological standard model. An enor- 
mous variety of observations at very large scales (> 1 
Mpc), from cosmic microwave background anisotropics, 
galaxy surveys, cluster abundances or Ly-a forest are 
successfully explained within this framework. However, 
to achieve these goals, the dark matter component is re- 
quired to behave as a collisionless non-relativistic particle 
fluid. This poses an important problem from the particle 
physics point of view. Indeed, no candidate with the as- 
sumed properties exists within the known particles, and 
therefore, physics beyond the Standard Model appears 
as the only feasible solution. Among the proposed can- 
didates, we find, on one hand, the axion which is the 
Goldstone boson associated to the spontaneous breaking 
of the Peccei-Quinn symmetry postulated to solve the 
strong CP problem of QCD. The production of axions 
in the early universe mainly takes place through the so 
called misalignment mechanism in which the angle is 
initially displaced from its equilibrium value = and 
oscillates coherently. Such oscillations can be intrepreted 
as a zero-momentum Bose-Einstein condensate which es- 
sentially behaves as a non-relativistic matter fluid. De- 
spite the fact that axions are light particles, this non- 
thermal mechanism produces cosmologically important 
energy densities. On the other hand we have the thermal 
relics, produced by the well-known freeze-out mechanism 
in an expanding universe. They are typically weakly in- 
teracting massive particles (WIMPs) such as the neu- 
tralino in supersymmetric theories |]| (for a recent re- 
view see also 0). In addition to their weak interac- 
tions with SM particles, these candidates usually also 
have very weak self- interactions (collisionless). 

Despite the success of CDM at large scales, the model 
exhibits certain difficulties at sub-galactic scales. In par- 
ticular, high resolution N-body simulations of dark halos 
show cuspy density profiles 0, which contradict ob- 
servations from low surface brightness galaxies |5j and 
dwarfs 0, Q that suggest flat density profiles. In addi- 



tion, CDM also predicts too many small subhalos within 
simulated larger systems [a, |2( , in contradiction with ob- 
servations of the number of satellite galaxies in the Lo- 
cal Group. Solutions which aim to reduce the power at 
small scales, but keeping the good properties of CDM at 
large scales, have been proposed at both the astrophysi- 
cal and a more fundamental level. They include modifi- 
cations of the primordial power spectrum at small scales 
lfj, or of the collapse histories of CDM and baryons 
ll| . but also modifications in the nature of the dark 
matter. There are two main proposals along the latter 
lines, namely, warm dark matter (WDM) |l2| and self- 
interacting dark matter (SDM) 13]. In the WDM sce- 
nario the particle dark matter is still collisionless, but 
much lighter with M < 1 keV. Free-streaming expo- 
nentially suppresses the power spectrum below the scale 
X FS ~ 0.2(fi Br /i 2 ) 1 /3(keV/M) 4 / 3 Mpc. This suppres- 
sion solves the sub-halo problem provided M ~ 1 keV 
|l4j . however the central cusp problem requires a much 
lighter particle M < 300 eV In addition, the Ly- 

a forest in quasar spectra observations puts a severe 
lower bound on the mass of the WDM candidate, namely, 
M > 750 eV, which makes the simplest models unfeasi- 
ble llflj l . A second possibility is to take into account dark 
matter self-interactions. Thus, if the mean free path of 
dark matter particles is of the order of ~ 1 kpc, but still 
keeping small annihilation cross sections, then heat can 
be transported out of the center of the halo smoothing 
out the central cusp, which now would become a more 
spherical core. In order for this mechanism to work, the 
total elastic dark matter cross section should be in the 
range (lif: 
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Ti 



— ~ 0.5 — 5 cm 2 



(1) 



which is comparable to a nucleon-nucleon interaction. 

However, the SDM model is also very strongly con- 
strained by different observations. Thus, in order to 
avoid evaporation of halos on a Hubble timescale, the 
cross section should be sufficiently small a/M < 0.3 — 1 
cm 2 g _1 01 • I n addition Chandra observations of the 
mass density profile of a certain cluster imposes a very 
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strong limit a/M < 0.1 cm 2 g^ 1 18]. Finally, the ellip- 
tical shape of a particular cluster determined by lensing 
data leads to a even stronger limit a/M < 0.02 cm 2 g _1 
[Tij. otherwise self-interactions would erase ellipticity. 
The interpretation of these results together with other 
direct constraints can be found in [gjj . It is also interest- 
ing to point out that the combination of both scenarios, 
i.e. self-interacting warm dark matter, could alleviate 
their individual problems [2l|. Apart from SDM and 
WDM, an interesting proposal in [22| shows that non- 
thermal WIMP production could also help resolving the 
mentioned discrepancies. 

The new requirements imposed by these models on the 
dark matter particles makes it even more difficult to find 
appropriate candidates in the mentioned SM extensions. 
For that reason it is worth analysing the potential of bra- 
non dark matter |23| as an alternative to the more com- 
mon neutralino or axion cases. In this paper we study 
the different possibilities in terms of the brane-world pa- 
rameters. We will consider the case in which branons and 
SM particles are the only relevant degrees of freedom at 
low energies. This happens whenever the brane tension 
scale / is much smaller than the fundamental scale of 
gravity in D = 4 + N dimensions Mp [24| . 

The paper is organized as follows: after introducing 
branon dark matter and its interactions, we study its col- 
lisional (collisionless) nature essentially in terms of Mn 
and /. Then we analyse the possibility of producing bra- 
nons non-thermally, in a similar way to the misalignment 
mechanism, where now the brane coherently oscillates 
along the extra dimensions. We end with some conclu- 
sions. 



BRANON INTERACTIONS 

In the case in which traslational invariance is an ex- 
act symmetry in the extra dimensions, the presence of 
the brane can break it spontaneously, and branons can 
be identified with the corresponding Goldstone bosons 
[25I Efil ]. For the sake of definiteness, we will con- 
sider a particular model with three compact extra di- 
mensions, although we expect that the results will de- 
pend only on the different energy scales involved in 
the problem and not on the details of the model. We 
will choose an extra space B = S 3 , so that its isom- 
etry group G{B) = SU(2) x SU(2) is spontaneously 
broken down to SU(2). Therefore, the coset space is 
K = SU(2)xSU(2)/SU(2) = SU{2) ~ S 3 = B. We thus 
have three branons, and the coset space metric is nothing 
but the 3-sphere metric: h a {3{ir) = S a p + ir a iT^ /(v 2 — ir 2 ), 
where v = f 2 R B with R B = M^ 1 (M P /M D ) 2 / N the size 
of the extra space B and Mp the usual Planck mass. 
Therefore, v is the spontaneous symmetry breaking scale 
which sets the size of the coset space K. 

As shown in |26j |. the low-energy branon dynamics is 



given, to lowest order in derivatives, by the non-linear 
sigma model corresponding to that symmetry breaking 
pattern, i.e.: 

C B r = ^"K^d^d^ + Q(p A ) 



+ 0(tt 6 ) + 0(p 4 ) 



(2) 



We see that v sets the strength of branon self-interactions 
to lowest order in derivatives. Notice however that the 
0(p 4 ) terms are not suppressed by the v scale, but by 
the / scale |26j. 

So far we have only considered massless branons, how- 
ever, traslational invariance in the extra dimensions can 
be explicitly broken, and branons can acquire mass. In 
such a case, we should include in the Lagrangian the cor- 
responding breaking terms. As shown in [27I Esl|. such 
terms will include an arbitrary number of even powers of 
branon fields, i.e.: 
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(3) 



Since branons are pseudoscalar particles, parity on the 
brane forbids terms with an odd number of fields |2fSj . 
Notice that here we are assuming for simplicity that the 
branons are degenerate, with a mass M which is small 
compared to v. Therefore the explicit breaking terms can 
be also organized in powers of the branon mass, which is 
considered to be of the same order as a derivative in the 
power counting. Accordingly, we expect the parameter 
A, which depends on the specific breaking mechanism, to 
be of order one. This is completely analo gous to the pion 
dynamics in Chiral Perturbation Theory 30] , in which v 
would play the role of the pion decay constant and the 
expression in J3J with A = 1 corresponds to the expansion 
in 7r fields of the lowest order mass term. 

Branons also interact with the SM particles through 
their energy momentum tensor. Again the lowest order 
term in derivatives was obtained in 26, 281 : 



\rpHV 

) A SM 



Cbt-SM = -^(Ad^dvTT** -M 2 Tr a TT a g^ + 



(4) 



where the dots stand for higher order terms in ir fields. 
We see that the branon-SM interactions are controlled 
by the brane tension scale /. What is important from 
the viewpoint of the dark matter problem is that bra- 
non self interactions are controlled by a completely dif- 
ferent energy scale from branon-SM interactions. This 
means that it is possible to have strong self-interactions 
and simultaneously small annihilation cross sections into 
SM particles, so that the branon relic abundance can be 
important. In addition, since branons always interact 
by pairs this implies that they are stable particles and 
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therefore they fulfill all the requirements of a dark mat- 
ter candidate. 



averaged over initial ones in the non-relativistic limit 
reads: 



COLLISIONAL VS. COLLISIONLESS BRANONS 

Branons were in thermal equilibrium with the rest of 
SM particles when the temperature of the universe was 
above their decoupling temperature, which depends on 
/ and M through the corresponding annihilation cross 
section. Assuming that the evolution of the universe was 
standard up to a temperature of the order of the brane 
tension scale, and that decoupling occured below such 
temperature, which is indeed the case (see 23]), then 
the above effective Lagrangians can be used to compute 
the branon thermal relic abundances generated by the 
standard freeze-out mechanism. The calculation of the 
abundances [23l |2jj show that there are two regimes in 
which the corresponding energy density could be com- 
patible with the measurements of the cosmological dark 
matter, i.e. flBrh 2 < 0.129 and also with LEP-II bounds 
pjj . On one hand we have a region in which M > 100 
GeV and / < M, which corresponds to cold relics, and on 
the other, a region with M < 180 eV and / > 200 GeV 
in which branons behave as hot or warm dark matter. 

In this section we explore these two regimes and calcu- 
late the branon self-interactions in terms of the number 
of branon fields N and the fundamental scale Md- With 
that purpose we will obtain that branon-branon elastic 
cross section. 

Since particles in the halo collide with very low relative 
velocities v re i ~ 10~ 3 , the scattering amplitudes can be 
obtained from the above effective Lagrangian and they 
are given, for A = 1, by the well-known Weinberg's low- 
energy theorem [3(| : 



Otot 



A(s, t, u) 



M 2 



(5) 



Assuming that the explicit breaking terms respect the 
remaining SU (2) symmetry, branons transform as an 
SU(2) triplet and it is possible to construct isospin 
I = 0, 1, 2 amplitudes as: 

To(s,t,u) — 3A(s, t, u) + A(t, s, u) + A(u, t, s) 

Ti(s,t,u) = A(t, s,u) — A(u,t, s) 

T 2 (s,t,u) = A(t,s,u) + A(u,t,s) (6) 

and also define the corresponding partial wave ampli- 
tudes with definite angular momentum J: 

i r 1 

tu = ^r d(cos6)P J (cosd)T I (s,t,u) (7) 

327T J_ t 

From the amplitude, the corresponding total elastic scat- 
tering cross section summed over final isospin states and 



23Af 2 
384ttw 4 



0(v 2 rel ) 



(8) 



Apparently this cross section grows unboundedly with 
M, this is due to the fact that we are working with an 
effective theory, and tree-level unitarity violations are ex- 
pected for sufficiently high values of s ~ AM 2 . In order 
to estimate the range of validity of the above result, we 
will use the well-known unitarity limit: 



Re tjj < 



1 



(9) 



where p = (1 — 4M 2 /s) 1 / 2 . From the dominant ioo wave, 
we get in the mentioned limit: 



M 2 16 7T 
— ^7 
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(10) 



which is compatible with the general result derived in 
31]. For typical values of v reU we find M/v < O(10 2 ). 
Substituting back in (JSJ, we find a limit for the cross 
section: 



< 4 ' 106 

M ~ M 3 



(11) 



In order for this kind of branons to be strongly self- 
interacting they should satisfy atot/M > 0.5 cm 2 g -1 = 
2.3- 10 3 GeV -3 , which implies M < 12 GeV, in agreement 
with the general bound 31]. This inmediately excludes 
cold branons as possible SDM candidates, according to 
the LEP-II limits mentioned before. Notice that, in any 
case, this is only a unitarity limit on the model, for larger 
masses, branons interactions could be sufficiently strong, 
but they will not be described by our effective theory. 

Nevertheless, self-interacting warm branons with 
masses in the region M < 180 eV are not a priori ex- 
cluded by the previous bound. In principle, such a light 
particle would conflict with Ly-a forest observations, as 
commented before. However, such limits were derived 
assuming collisionless WDM. In the collisional case free- 
streaming is suppressed and the relevant cut-off scale 
for structure formation is the standard Jeans length Aj, 
which can be substantially smaller than Xps UM- l n an y 
case a detailed calculation of the matter power spectrum 
would be needed for each particular model [2l|. 

In Fig. 1 we have plotted the different regions in the 
/ — Md plane for a branon mass M ~ 100 eV. The curves 
are not very sensitive to the branon mass, since they 
scale as M 1 / 8 . We see that for a given /, models become 
self- interacting at sufficiently high Md- Notice that the 
plotted regions are consistent with the assumption of low- 
tension brane / Md- 
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FIG. 1: Collisional vs. collisionless regions in / — Md space 
for a model with M = 100 eV and N=3,5,7. The dividing 
lines correspond to atot/M ~ 0.5 cm 2 g _1 . The lower striped 
area is excluded by LEP-II single photon analysis. 



THERMAL VS. NON-THERMAL BRANONS 

So far we have considered branons as WIMPs which 
were produced thermally in the early universe. However, 
if the evolution of the universe was not standard before 
branon decoupling or the reheat temperature Tpu was 
sufficiently low, branons were never in thermal equilib- 
rium with radiation. However, still they could be pro- 
duced non-thermally very much in the same way as ax- 
ions |3^ . Thus, if the maximum temperature reached 
in the universe was smaller than the freeze-out temper- 
ature Tjih -C Tf, but larger compared to the explicit 
symmetry breaking scale Trh ^> A with A = (Mi;) 1 / 2 , 
then brane fluctuations were initially essentially massless 
and decoupled from the rest of matter fields. In this case, 
there is no reason to expect that the position in the extra 
dimension Yq at which the brane is created should coin- 
cide with the minimum of the branon potential (Y = 0) . 
In general we expect Yn ~ O(Rb), i.e. ttq ~ v within 
a region of size H~ x |32|. The evolution of the branon 
field is then simply that of a scalar field in an expand- 
ing universe. Thus, while H 3> M, the field remains 
frozen in its initial position tt = ttq. Below the tempera- 
ture Tj for which 3H(Ti) ~ M, the branon field oscillates 
around the minimum. These oscillations correspond to 
a zero-momentum branon condensate, its energy density 
behaving like non-relativistic matter. 

For the sake of simplicity we will consider the case 
of branons which are lighter than all the SM particles 
masses, i.e. we will assume that M <C 1 MeV and that 
neutrinos are massless. In such a case the branon con- 
densate cannot decay by particle production into massive 
SM particles, but still branons could annihilate through 
2^2 processes into massless gauge bosons or neutri- 



nos. Since the expectation value of ir varies on horizon 
scales, the typical momentum of branons can be esti- 
mated as \p\ ~ H, i.e. \p\ ~ M at the beginning of 
oscillations (if there was a previous period of inflation 
then the typical momentum would be much smaller). It 
is possible to estimate the annihilation rate from the re- 
sults in [13 as Y 2 -, 2 ~ M 7 v 2 /f s at the first stages of 
oscillations. In the regions in which we will be inter- 
ested, we find T2^2 <C M. In addition, those processes 
in which four zero-momentum branons are transformed 
into two branons with energies E = 2M were already 
considered in the first reference in [32| . In this case 
Bose enhancement effects are important. The results, 
obtained for axions, can be inmediately translated into 
the branons case, simply identifying v with the Peccei- 
Quinn scale. It was shown that in the initial stages of 
oscillations, which corresponds to 7r ~ v, the decay rate 
was T^^n ~ M. Again, this rate would be smaller 
if there was a period of inflation. Therefore the nec- 
essary condition to avoid the energy depletion of the 
branon condensate is H(Trh) > M. Making use of 
the Friedmann equation in a radiation dominated uni- 
verse, this implies Trh > (MMp) 1 / 2 , which automat- 
ically ensures Trh 3> A. Since Y decreases very fast 
due to the universe expansion, then it will never be- 
come comparable to H. Accordingly, the branon con- 
densate energy density essentially is not reduced by par- 
ticle production, but only diluted by the Hubble expan- 
sion. We see that the reheating temperature should 
satisfy the condition T t ~ (MMp) 1 / 2 < T RH < Tf 
in order for the condensate to form and survive until 
present. Therefore if Tj > Tf the above interval dis- 
appears and only thermal relics can be present. In the 
opposite case, Tj < Tf, we can also have non-thermal 
production. In the case of light branons, a good esti- 
mation for the freeze-out temperature was obtained in 
23, 29]: log(7>/GeV) ~ (8/7) log(//GeV) - 3.2. In Fig. 
2 the Ti = Tf curve which separates the thermal and 
non-thermal regions is plotted as a dashed line. We see 
that for sufficiently light or weakly coupled branons the 
non-thermal production is possible. 

Following the same steps as in the axionic case, we can 
calculate the energy density which is stored today in the 
form of branon oscillations. Assuming that M does not 
depend on the temperature, it is given by [3^ : 

2 2.hNv 2 MT§ 6.5 • lCr 20 ^ 



17 Brh 



M P T p 



GeV 5/2 



-fR\M x l 2 



(12) 



where To and po are the photon temperature and critical 
density today respectively. 

In particular it can be seen that within certain param- 
eters regions, the above energy density can be cosmologi- 
cally important. Thus in Fig. 2 we plot the ng r h 2 = 0.1 
curves, corresponding to the current observational val- 
ues, for different values of Rb in the / — M plane. Thus 
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FIG. 2: Thermal vs. non-thermal branons regions in the / — 
M plane. The dashed (red) line separates the two regions 
and corresponds to Ti ~ T/. The continuous (black) lines 
correspond to ftsrh 2 ~ 0.1 for different values of the radius 



of the extra dimensions Rb = 10 , 1, 10 



10" 



GeV -1 . 



the region above each curve would be excluded by non- 
thermal branon overproduction. This is another man- 
ifestation of the cosmological moduli problem. Notice 
that in the allowed non-thermal regions of the figure the 
condition / <C Md is also satisfied. 

In the case in which the branon masses are larger than 
some SM particle mass, then the corresponding annihi- 
lation channel would open up and a detailed analysis of 
the energy loss rate would be necessary. In any case, the 
possibility of having heavier non-thermal branons is not 
a priori excluded. 



CONCLUSIONS 

We have explored the different possibilities of bra- 
non dark matter in terms of the brane-world parame- 
ters (N, Md, /, M). We have shown that apart from the 
standard scenario |23l | in which branons behave as colli- 
sionless thermal relics, the large parameter space of these 
models allows for branons to behave also as collisional or 
non-thermal dark matter candidates. In the collisional 
case, this opens an interesting possibility since the exist- 
ing proposals for SDM candidates in the literature are 
very limited [21L l33| . Nevertheless, the allowed parame- 
ters region reduces only to warm collisional branons and 
therefore a detailed analysis of the actual effect on struc- 
ture formation is needed. 
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